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We investigate the dependence of the Drell-Yan cross section on lepton polar and azimuthal angles, 
as generated by the lowest-order QCD annihilation and Compton processes. We focus in particular 
on the azimuthal- angular distributions, which are of the form cos (ji and cos 2(j}. At small transverse 
momentum qt of the lepton pair, qt ^ Q, with Q the pair mass, these terms are known to be 
suppressed relative to the </)-independent part of the Drell-Yan cross section by one or two powers of 
the transverse momentum. Nonetheless, as we show, like the ^-independent part they are subject 
to large logarithmic corrections, whose precise form however depends on the reference frame chosen. 

/— C ' These logarithmic contributions ultimately require resummation to all orders in the strong coupling. 

^—>. , We discuss the potential effects of resummation on the various angular terms in the cross section 

#vj ' and on the Lam- Tung relation. 
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I. INTRODUCTION 



(N 
> 

t~^ ■ It is well known that the angular distribution of the leptons in the Drell-Yan process P1P2 -^ £iX may possess 

azimuthal asymmetries. These correspond to angular-dependent terms in the ratio of differential cross sections, 

Oi dN _ r daY^ da 

^ ; dO = \d^) dndSi ' ^ -' 

where q is the four-momentum of the virtual photon (or Z boson, if the energy is sufficiently high) decaying into 
f^ • the lepton pair, and dQ = dcosOdcj) is the solid angle of the lepton £ in terms of its polar and azimuthal angles in 
I ' the center-of-mass system (c.m.s.) of the lepton pair. As we shall review below, an analysis of the general Lorentz 
M-i. structure of the hadronic tensor yields the following angular structure [ll 01 • 

^_ ■ dN _ S WTil+ cos^ e)+WL{l~ cos2 0) + Wa sin 261 cos </> + Waa sin^ 9 cos 2(j) 

'dn~8^ 2Wt + Wl ■ 

^^ , Here the "structure functions" Wt.l.a,aa depend on the virtual photon's invariant mass Q, its transverse mo- 
}^ ' mentum Qt, and its rapidity y. Within the lepton pair c.m.s. there is still freedom to choose the axes of the 
5t ! coordinate system, with respect to which the lepton angles are defined. The Wt,l, A, aa also depend on this 
choice. Equivalently to Eq. (0), one also often writes the lepton angular distribution as 



dN 3 1 



9 t/ 2 ' 

1 + A cos 6 + ^1 sin 26 cos (j) + — sin cos 20 



One obviously has 



dQ. 47r A + 3 

Wt-Wl Wa 2Waa 



Wt + Wl Wt + Wl Wt + Wl 



(3) 



(4) 
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A third, again equivalent, parameterization for the lepton angular distribution 3] that is also often employed is 

(5) 



dN _ 3 



l + cos2 6'+ —(l-3cos^ 61) +Aisin26'cos0+ — sin^ 6* cos 2(/) 



Evidently, 



or 



2-3Ao 2Ai 2A2 ,„, 

^ ^ 2Wl ^ ^ 2W^ ^ ^ m^^ 

° 2Wt^Wl ' ^ 2Wt + T^l ' ^ 2Wt^Wl ' 

If Q is large, one may calculate the structure functions using parton model concepts, factorizing them into 
collinear convolutions of the parton distributions of the two scattering hadrons and partonic hard-scattering cross 
sections that are amenable to QCD perturbation theory. For non-vanishing transverse momentum Qt of the 
virtual photon, the lowest-order (LO) partonic processes are qq — > 7*1? and oa— > 7*0. The contributions to the 
angular distributions by these processes have been calculated in Refs. 01 S H S 13 ■ Also the next-to- leading 
order (NLO) corrections have been derived |3ill3- 

Experimentally, for a given lepton pair invariant mass Q, the bulk of the Drell-Yan events is at rather low 
transverse momenta of the pair, Qt «C Q. It is this regime that is also most interesting from a theoretical point of 
view. "Intrinsic" transverse momenta of the initial partons may become relevant. Also, as Qt <C Q, fixed-order 
calculations of the partonic cross sections are bound to fail. When Qt -^ 0, gluon radiation is inhibited, so that 
only relatively soft gluons may be emitted into the final state. The cancellation of infra-red singularities between 
real and virtual diagrams in the perturbative series then leaves behind large logarithmic remainders of the form 
a^ In™ (Q^/Ql^) IQ\ in the cross section da/d'^q at the fcth order of perturbation theory, where m = 1, . . . , 2A:— 1. 
Ultimately, when Qt ^ Q, «« will not be useful anymore as the expansion parameter in the perturbative series 
since the logarithms will compensate for the smallness of as ■ Accordingly, in order to obtain a reliable estimate for 
the cross section, one has to sum up ( "resum" ) the large logarithmic contributions to all orders in as ■ Techniques 
for this resummation are well established, star ting with pioneering work mostly on the Drell-Yan process in the 
late 1970's to mid 1980's [H IH El El El Ei IH Eg. The "Collins-Soper-Sterman" (CSS) formalism [3 
has become the standard method for Qt resummation. It is formulated in impact-parameter (6) space, which 
guarantees conservation of the transverse momenta of the emitted soft gluons. 

The large terms a^' In™ {fP /Qt) IQ\ just described occur only in the part proportional to Wt of the cross 
section da/dfld'^q. The function Wa is less singular by one power of Qt at small Qt, while Wl and Waa are 
suppressed even by a factor Q^ relative to Wt- For Wa and Waa this can be understood as follows: at Qt = 
the angle (f> cannot be defined, hence no azimuthal asymmetry could be observable, and therefore, the asymmetries 
should smoothly go to zero in the limit of Qt -^ 0. 

If one is just interested in the small-Qr behavior of the cross section da/d q, it will be sufficient to take into 
account the resummation of the large logarithms in Wt- However, this may become different if one considers the 
parameters A, /i, v or ^0,1,2 defined above, which were the object of dedicated experimental studies in Tr^iV Drell- 
Yan experiments about 20 years ago El |2l ^^ (experimental evidence for nonzero Wl was already reported 
in [2j|). Superficially, one might think that, in order to improve the theoretical prediction for these coefficients 
by resummation, it will be sufficient to perform the resummation in Wt alone. This is in fact what has been 
done in the literature so far in |23| for Drell-Yan and in [2J, l25L .26.] for the related "semi- inclusive deeply- inelastic 
scattering" (SIDIS) process and has led, for example, to the claim in Ref. [23 that resummation has a strong 
effect on the pcrturbative-QCD results for A,/^ and v at small Qt- However, even though Wl, Wa, and Waa 
are down by powers of Qt as Qt —^ 0, they all individually may receive very similar large logarithmic corrections 
at small Qt as Wt does. Therefore, it may well happen that in a ratio such as A2 = 4Waa/(2Wt + Wl) the 
resummation effects cancel to a large degree, if not completely. In the present paper, we will investigate the 
small-Qr behavior of Wl, Wa, and Waa, and we will also address some qualitative consequences for the Qt- 
resummation for these and for the various coefficients constructed from them. A complication arises from the fact 
that the Wl, Wa, and Waa depend on the coordinate frame chosen. In a change of frame, terms proportional 



to Qt or Qy may be redistributed among Wt and the M^l,a,aa and, because Wt is more singular at small Qt, 
may alter the small-QT behavior of Wi_A,AA- Our aim is to exhibit this feature very explicitly by considering 
two particular frames commonly used in the literature, in order to eliminate some misconceptions concerning the 
small-QT behavior and the resummation of angular distributions. 

As discussed above, the CSS formalism [T^ describes the behavior of the structure function Wt at small Qt- 
It may be used to predict the Qr-singular pieces arising at a given perturbative order. To LO, these have a 
particularly simple structure [13, involving the first-order expansion of the Sudakov form factor and the LO 
DGLAP ^3 splitting functions. This structure is straightforwardly recovered from explicit calculations of the 
cross sections for the partonic reactions qq — > 7*3 and qg — > 7*g. When we use these processes to determine the 
small-Qr behavior of the other structure functions Wl, Wa, and M^aa, we find a closely related, but different, 
form that involves different "splitting functions" . This is likely due to the fact that the direction of the observed 
transverse momentum matters for these structure functions. It is important to emphasize that the CSS formalism 
was not constructed to treat the directional dependence of the transverse momentum distribution. Even though 
we will discuss in this paper some features of the resummation of the Wl^a.aa, we have not been able to organize 
their small-Qr behavior beyond the leading logarithms into a resummed form in terms of an "extended" CSS 
formalism that goes beyond coUinear factorization. We therefore hope that the study presented in this paper will 
serve as a motivation for a more general analysis of resummation in cases where the direction of the observed 
transverse momentum matters. In this respect the proper starting point will be to consider the factorizations and 
techniques put forward in Refs. [l5llla.l2q. 

The article is organized as follows. In Section II we express the azimuthal asymmetries in terms of structure 
functions in two different frames that are commonly used in the literature. This analysis is exact in Qt and does 
not rely on the use of perturbative QCD, just on kinematical considerations. In Section III we discuss the LO 
perturbative QCD results, and in Section IV we investigate their small-Qr limit. In Section V we address some 
aspects of the resummation of the angular terms in this limit. We present our conclusions in Section VI. 

II. AZIMUTHAL DEPENDENCES IN TERMS OF STRUCTURE FUNCTIONS 

In this section we discuss in detail the (p dependence of the differential cross section in terms of the structure 
functions Wt.l, A, AA- This is a purely kinematical analysis, which will set the notation to be used later on. We 

follow (in part) the notation of Refs. 2e2M^ tiut use the metric diag(H ). In passing we will point out some 

differences with results that appeared in the earlier literature. 

We first define invariant structure functions 1^1,2,3,4 through the following parameterization of the hadronic 
tensor: 

q'' 2 



y^^u ^ _^g^u _ 1^)W^ + P^^P''W2 ~ -{P^p" +p^'P'')W3+p^'p''W4 , (8) 



where P = P^ + P^, p = Pi - P2, P'' = (P^ - q^'{P ■ q)/Q'^)/Vs, p" = [p" - q^ip ■ q)/Q^)/^Vs, with Pi,2 the 
initial hadron momenta and q the momentum of the virtual photon, Q^ = q^. The hadronic c.m.s. energy is 
Vi=v/(A+^2)2. 

As mentioned in the Introduction, we will use two different reference frames: the so-called Collins-Soper (CS) 
frame Q and the Gottfried- Jackson (GJ) frame J4I (also sometimes referred to as the "i-channel helicity frame"). 
Both frames are rest frames of the virtual photon or, equivalently, of the lepton pair. However, this property 
does not completely specify a frame, as one has a freedom in the choice of the coordinate axes, corresponding 
to rotations among the frames. In both the CS and GJ frames, a set of orthonormal axes {X, Y, Z and T) are 
defined, and the hadronic tensor is re-expressed as: 

W"" = -{g"" - T''T''){Wt + Waa) - 2X''X''Waa + Z'^Z'^iWi -Wt- Waa) - {X^Z" + Z''X'')Wa ■ (9) 

In this way one has W^f^ = —{2Wt + Wl), which is the quantity that appears in the differential cross section 
da/d'^q, the denominator of dN/dil in Eq. |^. Contracting with the straightforwardly calculated leptonic tensor, 
one finds the cross section in terms of the structure functions: 



da c? 



dnd'^q 2(27r)4g2s2 



{WT{l + cos^e) + WL{l~-cos^e) + WASm2ecos(j) + WAAsin^ecos2(l)} , (10) 



again valid in both frames. Here a is the electromagnetic coupling constant. From Eq. H10() one immediately 
reproduces Eq. ^. Higher harmonics in cos{n(j)) do not occur in the angular distribution due to the fact that 
the kinematics of the process are fully determined by only three momentum vectors, Pi, P2, and q. In case the 
polarization of the leptons is not summed over or in case of electroweak corrections sin cj) and sin 2(j) terms will 
also be present, but these will not be considered here. We also note that the structure functions Wt, l, A, aa are 
associated with specific polarizations of the virtual photon 0: Wt = W^^'\ Wl = H^°'°, Wa = {W°''^ + W'^'°)/V2, 
and Waa = W^'"^, where the first (second) superscript denotes the photon helicity in the amplitude (its complex 
conjugate) of the Drell-Yan process. The azimuthal dependence introduced by Wa and Waa therefore comes 
from their single- or double-spin-flip property, respectively. 

In both frames, T^ = g^/Q and V^ = e^^'^^^Xi^ZaTp. In the CS frame, the Z axis is defined as pointing in the 
direction that bisects the angle between the three- vectors P2 and —Pi, see Fig.^ This gives 



Z^ = 



VQ^+Q 



(qpP^ + qpp^) = 



*V^ 



^T 



((P2 • g)Pr - (Pi • q)Pi 



(11) 



where Q^ = q^ is the square of the transverse momentum q^ of the virtual photon with respect to the two 
hadron momenta. Furthermore, qp = (P ■ q)/y^, qp = —(p ■ q)/\/s (in the hadronic c.m.s. qp is the energy of the 
virtual photon, while qp is its ^-component). Finally, in Eq. (|11() P^ = P^ — q^{Pi ■ q)/Q^ (note there is no factor 
1/y/s in the latter definition, compared to the definitions of P and p above). For the X-axis one chooses 



X'' 



Q 



^VQ' 



{qpP'' + qpp^) 



2Q 



sQtVQ^ 



Qrp 



((P2-g)Pr + (Pi-g)Pr). (12) 



In the GJ frame, the Z axis points in the direction of the three- vector Pi, see Fig.|21 



Z^' 



qp -qp 



P^" +p^ 



^p. 



Pl-q 



(13) 



Furthermore, 

X^' -- 



2Q 



QriQ^ + Qj^) 



{{Q^qp - Qlqp)P^ + {Q\p - Qlqp)p^) = -^ ((P2 • Z)P!^ - (Pi • Z)P.^) .(14) 



One finds that the three-vector components of the X and Z axes of the two frames are related by a rotation: 



Zgj = cos 7 Zcs + sin 7 Xcs ■. 
Xgj = - sin 7 Zcs + cos 7 Xcs 



(15) 
(16) 




lepton plane (cm) 



FIG. 1: The Collins-Soper frame. 




lepton plane (cm) 



FIG. 2: The Gottfried- Jackson frame. 



where 



cos 7 ■ 



Q 



VO^ 



fr 



sin7 



Qi 



^Q^ + Ql 



(17) 



Thus, in the Hmit Qt — > the two sets of X, Z axes coincide. 

Inserting the above definitions of the coordinate axes into the hadronic tensor in Eq. (Q and comparing to 
Eq. ((21), one can derive for each frame the relations between the sets Wi,2,3,4 and Wt,l, A, aa of structure functions. 
For the CS frame this gives (this is Eq. (B2) of Ref. 0): 



Wt = Wi- l^AA 
1 



Wl = Wi 
Wa = - 
W^AA - - 



Qt 



{qpW2 + qpqpWs + qlWi) 
1 



g(Q2 ^ g2y \ qpqp{W2 + Wi) + ^{ql + ql)W3 ) , 



Q', 



2Q2(Q2 + Qi,) 



{qpW2 + QpqpWs + q^Wi) 



(18) 



while for the GJ frame one finds (this is a corrected version of Eq. (B3) of ^ and also of a corresponding expression 
in reference-note 5 of |2Sy): 

Wt = Wi~ Waa , 
Wl = Wi + {a^W2 + af^Ws + (3^Wi) , 
Qt 



Wa 



qp - qp 



Waa = - 



Qt 



2('?p - qp)'- 



a{W2 + Iws) + (3{W4 + I-W3) 



{W2 + 1^3 + W4) 



(19) 



where 



a 



qp 



Q 



Q qp - qp 



P^'^- 



Q 



Q qp - qp 



(20) 



The relation between the structure functions in the CS frame and the GJ one is given by the foUowing hnear 
transformation: 




p2 



/I 

V 



p 


\/ 


-p 


\/ \ 




/ Wt 


1 


2p 


-p" 




Wl 


p 


-p 


l-p' 


-P 




M^A 


¥' 


-¥' 


P 


1 + ipV 




Vm^aa 



(21) 



GJ \ -If" 2^ F - i 2^^ / x".^"/ CS 




where p — Qt/Q- Expressed in terms of the angle 7 defined above, the matrix in the above equation is identical 
to the one presented in |29| (their Eq. (4)), apart from a sign in the third entry of the first row. 

For future reference, we also give the transformation between the sets of coefficients A,/i, v in the two frames, 
which were defined in Eq. Q). From Eq. (|21() one finds: 

-ip |p2 \ /A\ 

(22) 

GJ \ I" ">" ^ ' If" / X' / CS 

where 

A = l + p2 + ip2A + p^-ip2^. (23) 

The reverse transformation from the GJ frame to the CS frame is the same upon replacement of p ^ —p (and 
exchange of the labels CS and GJ^. This is in agreement with Ref. J,9] (note that what is referred to as "m- 
channel" (UC) frame in that reference corresponds to the GJ frame used here, which accounts for the p — > —p 
difference with respect to our Eq. H22(l ). One notes that if p — > 0, the rotation matrix becomes the unit matrix, 
as expected. All results presented so far are exact in p. 

III. AZIMUTHAL ASYMMETRIES FROM LO PERTURBATIVE QCD PROCESSES 

The structure functions can be calculated employing coUinear factorization and QCD perturbation theory. One 
writes down an expression analogous to (|HJ or |(^ for the partonic tensor, in terms of partonic structure functions 
Wi. The hadronic structure functions are obtained as convolutions of the partonic ones with the appropriate 
parton distribution functions. The Wi may be evaluated directly from partonic hard-scattering processes, which 
at LO are the annihilation (or, gluon bremsstrahlung) reaction qq — > 7*3 and the (QCD) Compton process 
qg —f 7* (7. The annihilation process is the dominant process in n^p Drell-Yan scattering (likewise in the pp 
Drell-Yan process). The Compton process is dominant in the Drell-Yan process in pp scattering at large Qt- 

In the following, we will first focus on the annihilation process, which for our purposes is most relevant, since the 
Compton process is sub-leading in the region of Qt <C Q we are mostly interested in. For the process qq — > 7*17 
one has: 

a;iX2 J J ' ^ i\t,2 

where we define X\ and Xi by writing q = X\P^ -I- X2P^ + qr, with the light-cone components of any four-vector 
V given by v^ = {v^ ± v^)/\/2. We have xi = {qp + qp)/^/s, X2 = {qp — qp)/y/s with qp,qp as introduced in 
the previous section (note that in the hadronic c.m.s. one has qp = {s + (3^)/2^/s, qp = \/ q^p — Q^ — Q\ — 
\/{s — Q'^Y I'^s — Q^ ). For the partonic coUinear momenta we define pf = £,iP^ and p^ = £,2P2 ■ The delta- 
function in if^ expresses the on-mass-shell condition for the outgoing "unobserved" gluon in the process qq -^ 7*5. 
Equation (Pljl contains the appropriate sum over all quark and anti-quark flavors a, each with their corresponding 
parton distributions qaix, p) and qa{x, p) = qa(x, p), respectively, where p ^ Q is the factorization scale; e^ is the 
quark's squared electromagnetic charge. Note that since we are only interested in the lepton angular distribution 



dN/dfl, which is a ratio of cross sections, we are free to adjust the overaU norniahzation of the tensors, which we 
do in such_a way as to simphfy the formulas. The annihilation process then has the following partonic structure 
functions 



W^ = 



W2 



W. = 



Cfs 



2tt Ci^Qt 
as CfQ^ 

Cpas Q^ 
2 

T 



Q 2Qt t- f 



efCl 



[e.^^i] 



as Cf sQ"^ 



27r til 



ZTT tu 



U2Q 



n~&i-'^''-^^n~e.] 



2-n 



tu 



Wi 



Wo 



(25) 



where Cp = 4/3, and where s ~ (Ci A + ^2-^2)^ = C1C2S, i = (q — £.iPi)^ and u — {q — £.2P2)^ ■ Constructing via 



Eq. I|18() the structure functions Wt, l, A. AA and inserting these into Eq. (0), one finds 



in the CS frame: 



dN 

'dn 



3 

16^ 



9 ^X' ^ 9 ^T 



COS 9cs 



QtQ 
Q2 



■K{xi,X2, Qt/s) sin 20cs cos ( 



1 



'CS 



?T 



2 Q2 + Q 



sin 9csCOs2(l)cs 



where the function K{xi, xi-, Qt/s) is given by 

/d6 /rf6 5 ((a - Xi)(6 - X2) - Ql/S) Eae^'Za(6,M)ga(6,M) H/^? - ^^l/Cl) 



X(a;i,X2, Qt/s) 



(26) 



(27) 



/d6/d6'5((a-a;i)(6-2:2)-Q|/s)Eae^'?a(a,M)9a(6,M)(a;?/ef+a;l/e|) 
As one can see, for the annihilation contribution in the CS frame all effects of the partonic light-cone momentum 
fractions and the parton densities cancel in dN/dfl, except for the term involving sin 29cs cos (j)cs associated with 
the ratio Wa/{2Wt + Wl)- As a consequence, the coefBcients A and ly defined in Eq. ^ are also free of any 
dependence on the parton distributions to this order Q. One reads off: 



Acs 



Q' ~ hQ 



Q2 + m 



vcs 



Tt 



Q^ 



^Qt 



(28) 



so that 



1 - Acs - 2j^cs = , (29) 

which is the well-known Lam- Tung (LT) relation |2,l23- I^ i^ equivalent to Wl — 2Waa and Aq = A2: its origin 
is the equality of W2 and W4 in Eq. f^ . 

In the Gottfried- Jackson frame one finds that all terms in dN/dVt depend on the parton densities and do not 
have particularly simple expressions. Remarkably, the LT relation continues to hold in this frame, however. This 
may be seen from Eq. H22(l which readily shows that the LT relation indeed holds in the GJ frame if it holds in 
the CS frame, and vice versa. In fact, it turns out that the LT relation holds for any definition of the lepton pair 
c.m.s. frame to this order. 

7*(7 subprocess. The contributions to the partonic structure functions for 9(^1)5(^2) 



Next, we consider the qg 
are Q: 



Wl = 



W2 



w. 



Wi = 



as 2Tr 

2-K si 
as 2TrsQ^ 
27r si 
as ATrsQ^ 
2tt si 
as 2TrsQ^ 
27r si 



[{S-Qr + {i-Q')'] 

[26(^1 + 6) + el] , 

[el-2e?] , 

[26(^1 - 6) +^2] , 



(30) 



where Tji — 1/2. For the contribution from .9(^1)5(^2) one has to interchange ^1 and ^2, i and u, and, in addition, 
change the sign of W3. 

From the expressions in H3()(l one finds that for the Compton process the ratios A, fi, v all depend on the 
parton distribution functions. As discussed in Ref. [Sj, there is in fact no frame in which any of the azimuthal 
distributions for this process become independent of the parton densities (only for very large Qt does this become 
approximately the case in the CS frame). Even if there were a frame in which this happened, the result would 
have no practical relevance, since one would always need to add the contribution by the annihilation process in 
the numerator and the denominator of the azimuthal ratios, which would spoil the cancellation of the parton 
distributions anyway. In fact, even the cancellation of the parton densities in Acs and vcs in Eq- H28|) for the 
annihilation process is of limited use, except in a pure flavor non-singlet situation where the Compton process is 
absent. ^^ ^^ 

However, despite the fact that W2 7^ W4, the LT relation does hold for the Compton process as well, and 
therefore for the complete LO Drell-Yan cross section \3L l30l| . In the next section we will verify the LT relation 
very explicitly in the small-QT limit. However, it is valid at LO regardless of the value of Qt- It is known to 
be (mildly) broken by NLO corrections [i| [l(J . Experimentallv the LT relation was found to be rather strongly 
violated in the 7r~p Drell-Yan process |l^ [23, UM ' ^ disagreement (both in ma gnit ude and in sign ) w ith the slight 
violation predicted at NLO. This has prompted much theoretical work |3lL 133. ISa. l34L ISa . ISa. 1371 IsM l39 | . offering 
explanations that go beyond the framework of coUinear factorization and perturbative QCD to which we restrict 
ourselves in this paper. 

As a final point in this discussion of the LO contributions, we would like to stress that the above discussion is 
not specific to the Drell-Yan process, but also appHes to similar processes like SIDIS |2J, [23, lla, 113 and back-to- 
back hadron production in two-jet events in electron-positron annihilation 41], when appropriate definitions of 
the frames and coordinate axes are used. 

IV. small-Qt limit of the structure functions Wt,Wl,Wa, Waa at LO 

The LO results presented in the previous section should provide good approximations for the angular distribu- 
tions in the Drell-Yan process at large transverse momentum, Qt ^ Q- We shall now investigate their behavior 
at small Qt- This is extracted most conveniently by using the expansion |24| 

6 ((1 - zi)(l - Z2) - Q't/§) = ^r—^ + 7T—^ '^(l - ^i)'^(l - ^2) lnp2 + 0{p') (31) 

for the delta-function in Eq. H24|) . where Zi = Xi/^i. Here, the "plus" -distributions are defined as usual for an 
integral from a; to 1 as 



2 (1 — z)_|_ Jx z 1 — z 

for any suitably regular function /. Inserting the structure functions Wi of Eqs. (|25|) and (|30|) into Eqs. (|18|) . and 
expanding for small Qt with the help of l|31|l . one finds in the CS frame: 

Wt,cs = ^— [-Cf (21np2 + 3) q{xi)q{x2) + q{xi) {Pgg ® q) (xa) + {Pqg (E) q) {xi)q{x2) 



2tt p 



Wl,cs = 2Waa,cs , 
a 



+ q{xi) {Pqg g) (a;2) + {Pgg g) ixi)q{x2) + ©(p^)] , 



2^ [-^F (21^/0^ + 3) q{xi)q{x2) + q{xi) {Pgg ® q) {X2) + {Pgg ® q) {xi)q{x2) 



Wa. 



CS 



as 1 
27rp 



+ q{x,) {P'^g ® g) {X2) + {P'gg ® g) {xi)q{x2) + 0{p')\ 
q{xi) [Pgg (g) qj {x2) - [Pgg ® q) {xi)q{x2) 



+ q{xi)(^Pqg^g'^{x2)-[Pqg^9){xi)q{x2)+0{p^)'\ , (33) 

where in each case we have kept the terms with the leading-power behavior at small Qt- Furthermore, we have 
for notational simplicity suppressed the factorization scale /z ~ Q in the parton distributions, as well as the sums 
over flavors, and we have defined the usual convolutions 

ir<S>f){x,)^ f'^V{x)f(^) , (34) 

with V variously one of the well-known ^7\ LO splitting functions 

' 1 + x^ 3 



P,,ix) = C_ 
or one of 



(l-x)+ 2 



:^'^(1-^) 



, P,g{x)^TR[x^ + {l-x)^] , (35) 



P:,gix) = P,g{-x) , Pgg{x) = Cf{1 + x) , Pgg{x) = Tfl ( 1 - 2x^) . (36) 

Note that one could cancel the term —3CFq{xi)q{x2) in Wt.cs against the contributions by the 3CfS{1 — x)/2 
term in the splitting function Pqq in Eq. H35(l . We have however kept the term in order to have the full splitting 
function in H35|) . and also because the term — Ci? (2 In p^ -|- 3) is the well-known first-order contribution to the 
Sudakov form factor. We also note that we could have first taken the small-Qr limit of the structure functions 
^^1,2,3,4 in Eqs. 123), l|^ and then inserted the result into (fT^ . In that case, since all the Wi have the overall 
power Qy and may be multiplied by powers of Qt in the transformation H18|l , it would have been crucial to keep 
the first sub-leading term proportional to p° in the Wi. Otherwise, one would obtain an incorrect result. 

As can be seen in Eqs. H33(l , the annihilation process makes a logarithmic contribution to each of the structure 
functions Wt, Wl, and Waa, on top of their nominal power in Qt- The Compton process, on the other hand, 
does not produce this leading logarithmic behavior. It is also interesting to note that, unlike the other structure 
functions, Wa,cs does not receive a logarithmic contribution at all. This turns out to be a result specific to the 
CS-frame. In case of Wt, the structure at small Qt is well- understood in terms of the CSS formahsm [T^l, as we 
shall briefly review in the next section. The LO small-Qr expressions for the other structure functions are new. 
We note that the Lam- Tung relation Wl — 2Waa of course still holds in (|33|l . 

The small-QT expressions for the structure functions Wt. l. A, aa in the GJ frame can be obtained in the same 
way, by using Eqs. 1)19(1 . or alternatively Eqs. 1(21(1 . and we find: 

Wt,gj — Wt,cs , 
Wl,gj = 2Waa,gj 

= 2^ [~Cf (21np2 + 3) q{xi)q{x2) + q{xi) (P+ ® q) (xs) + (P,- ® q) {xi)q{x2) 

2\ 



a 1 

Wa,G.J = TT- h^^ (21np^ + 3) q{xi)q{x2) + q{xi) {P+ (g) q) (xs) + {Pgg ® q) {xi)q{x2) 



271 p 



where 



+ 



P±(X) = Pgg{x) ± Pgq{x) 



+ 9(2^1) {P!,^ ® g) (^2) + {P;,- ® g) ixi)qix2) + Oip')] 
>+ (g> q) (X2) + {Pg- 

qixi) (P+ ® g) (xs) + [p-g ® g) {xi)q{x2) + 0{p^)\ , (37) 



qq y-^J — ^ qq\-^J -^ -^ qq\ 

P't{x) ^ Pgg{x)+P'{x)±2Pqg{x) 



(2Tr{1-x) 

P^gix) ^ Pqg{x)±Pqg{x)^\ (38) 

I 2ijfX(2x — 1) 
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with Pqq{x), Pqq{x), Pqg(x), Pqg(x), P' Jx) as given above. As one can see, in the GJ frame, all structure functions 
receive logarithmic contributions at small Qt- We observe that apart from Wt,gj none of the functions contains 
sub- leading terms (i.e., terms non- logarithmic in p) that involve only the usual splitting functions Pqq and Pqg of 
Eq. l(55|l . Therefore, Wl and Waa are not proportional to Wt at sub- leading order for small p, not even if one 
restricts to the annihilation process. The LT relation holds as before. 

Comparing the small-Qx behavior in the CS and GJ frames we find that the rotation between the two frames 
simply "reshuffles" the various "splitting functions" appearing in the Wt_l^a,aa- For example, in the annihilation 
contribution the functions contributing to the sub-leading terms in Wl_a,aa are Pqq and ±Pgq in the CS frame, 
but Pqq ± Pqq in the GJ one. This pattern, which also extends to the Compton part, can be verified by inspection 
of Eq. laj. 

We finally note that the non-standard "splitting functions" in the above expressions for the small-Qr limit 
are associated with the polarization states of the virtual photon contributing to the various structure functions. 
Only if the photons in the amplitude and its complex conjugate both have the same, and physical (transverse), 
polarizations, does one recover the ordinary DGLAP splitting functions at small Qt- This is only the case for 
Wt, cf. Section II. 

V. EFFECTS OF RESUMMATION AT SMALL Qt 

The small-Qr behavior of Wt is predicted by the CSS formalism ^M, which resums its singular behavior at 
Qt IQ ^ to all orders in the strong coupling constant. The CSS resummation is formulated in impact-parameter 
space. Schematically, keeping only those terms in the formalism that play a role for our present study, one has 

W^T = I ^ '''"' T. ^' 9a(xi, 5o/6) Qa(x2, 60/6) e^^^'^^) . (39) 



Here, S'(6, Q) is the Sudakov form factor, given by 

Ibl/b^ ^T 



S{b,Q) = -l_ ^ 



9" 

kl 



A{as{kT))\n(^]+B{as{kT)) 



(40) 



where 60 = 2e '^'^ with je the Euler constant, and where the functions A and B have perturbative expansions of 
the form 

A(a.) = 5]A.(^) , B{a,) = Y^B,[^)\ (41) 

fe=i fc=i 

Here we step over the fact that the b integration also may require the introduction of a nonperturbative Sudakov 
form factor, and simply focus on the perturbative part. 

For our purposes we only need the coefficients Ai and Bi, which read: 

Ai^Cf, Bi= -^Cf ■ (42) 

The term oc Ai in the Sudakov exponent generates the leading logarithms, which in &-space are of the form 
ajln^ (bQ), corresponding to a^ In^ ~^ {Q^/Qt) /Qt ^^ Qr-space. Next-to-leading logarithms (NLL) are gen- 
erated by the term oc i?i, and by the running of the strong coupling and of the parton distribution functions 
in H39() . In the CSS-resummed expression, the factorization scale is bo/b. Using the customary DGLAP evolution 
equation, the parton densities at this scale may be expressed by their values at scale fi ^ Q. Expanding the 
corresponding evolution matrix for the parton densities and the Sudakov exponential in as and considering only 
the first-order term, we recover the expression for Wt,cs given in Eq. H33|) (or \67]i ). The terms of Wt.cs in H33|) 
that involve the gluon distribution are generated by the singlet mixing in the evolution of the parton distributions 
between scales bo/b and Q. Clearly, for Wt,cs in Eq. 133|l to be reproduced by the CSS formalism, it is crucial 
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that the sphtting functions in its expression are the usual DGLAP sphtting functions [23 Pqq and Pqg that are 
associated with the evolution of quark and anti-quark distributions. 

This implies that the resummation of the structure functions Wl, Wa, and Waa will not precisely follow 
the CSS formalism. Taking their expressions 1331) in the CS frame as an example, we notice that the first-order 
expansion of the Sudakov form factor appears in Wl^cs and Waa,cs, indicating that the rcsummed expressions 
for these will contain the exponential of S{b,Q), as the one for Wt does, and hence have the same leading 
logarithms. However, the "splitting functions" involved in the non-logarithmic pieces are not the usual ones, as 
we already observed in the previous section. For Wa,cSi even the Sudakov part is absent, and the "splitting 
functions" are different again. All this means that beyond NLL the resummed expressions for Wl,c'S, Wa,cs, 
and Waa.cs will not organize into the structure given in Eq. H39(l . They likely will have a similar structure, but 
contain additional terms depending on the direction of b. To be clear, this is not a problem: while the presence of 
the true spitting functions in the leading term Wt at small Qt is required by (coUinear) factorization, this is not 
the case for the other structure functions. In other words, the "splitting functions" in these are not associated 
with standard parton evolution or coUinear singularities. Therefore, the structure we find only means that the 
resummation of the Wl.cs, Wa,cs, and Waa,cs (and likewise in the GJ frame) is more complicated beyond the 
leading logarithms. We do not address their full NLL resummation in this work, but note that the techniques of 
Refs. UM 13123 th^-t SO beyond coUinear factorization should prove useful for this purpose. 

Even without performing the full resummation, we can make some qualitative observations regarding the effects 
of resummation. The first observation from Eqs. H33I) and (|37|l is that the LT relation will not be affected by 
resummation. Next, we focus on the angular coefficients A,/i,t^ or Aq, Ai, A2 introduced in Eqs. (O and lO, 
respectively, which are ratios of the structure functions. Because of the structure of the leading logarithms, we 
expect that in the CS frame A and i/, and ^0 and A2, will be rather insensitive to resummation effects. Note that if 
the annihilation process alone contributes (which could in principle be realized by considering a flavor non-singlet 
combination of cross sections), resummation effects cancel identically in Acs and vcSi and they retain their LO 
forms given in Eq. (j28|l even after resummation. On the other hand, /i and Ai, for which there are no leading 
logarithms in the numerator because of their absence in Wa,cSj will be subject to substantial modification by 
resummation. In the GJ frame all ratios A, n, v will be affected by NLL resummation effects. These need not be 
small, in particular if the overall effects of resummation on Wt are themselves large. 

Next, we confront our findings with results of the previous literature. Chiappetta and Le Bellac [23 were the 
first to study the effects of resummation on the azimuthal asymmetries in the Drell-Yan process. Working in the 
CS frame, they argued that since the structure functions Wl,Wa and Waa are less singular than Wt, only the 
latter requires resummation. Therefore, they took into account resummation only in the denominators of the Ai, 
while for the numerators they employed the LO expressions. In this way they found large effects of resummation 
on the Ai. In light of our discussion above, the neglect of resummation in the numerators of the Ai is not justified. 

Similarly, in studies |2j, |2^ |2^ of resummation effects in semi-inclusive DIS only the resummation of the 
^-independent terms was taken into account. Specifically, resummation was not applied in studies of the ratio 
(cos2(/))/(cos(/)). Here, a Gottfried- Jackson type of frame (one of the hadron momenta defined the Z-axis) was 
used. Our results above show that the numerator and denominator of this quantity will have a somewhat different 
resummation, even though it may well turn out to be the case that the residual effects on the ratio are small. It 
would be desirable to revisit this quantity in the framework of a full NLL resummation study. This also applies 
to azimuthal asymmetries in polarized scattering J42|. 

VI. SUMMARY 

We conclude by summarizing our main observations. We have studied the structure functions Wt. l. A. aa in 
the Drell-Yan process, on the basis of the contributions by the annihilation process qq — > 7*17 and the Compton 
channel qg — > 7*g. The structure functions Wa and Waa in particular generate azimuthal asymmetries in 
the angular distribution of the produced leptons. The structure functions depend in general on the choice of 
coordinate axes, and we have investigated the results in the ColUns-Soper and the Gottfried- Jackson frames. 

We have focused on the behavior of the structure functions at small transverse momentum Qt of the lepton 
pair. We have recovered the known small-Qj- behavior of Wt, in terms of the leading-order DGLAP splitting 
functions and of the first-order expansion of the Sudakov form factor. For the other structure functions, which 
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are all nominally suppressed by one or two powers of Qt/Q with respect to Wt, we have found that they, too, 
in general have large logarithmic terms at small Qt, whose form depends on the frame chosen. We are not 
aware that this feature was pointed out previously in the literature. The small-Qr structure we find for Wl^a.aa 
differs from that of Wt, however. In the CS frame, Wl and Waa receive large leading-logarithmic corrections 
identical to the ones in Wt- By contrast, these are absent in W^. In the GJ frame, all structure functions have 
leading-logarithmic terms. Both frames have in common that the subleading terms in Wl^a,aa are different from 
those in Wt- This implies that the next-to- leading logarithmic resummation at small Qt must proceed differently 
from that for Wt given by the CSS formalism which uses a coUinear expansion. 

Without actually deriving the NLL resummation, we have discussed some generic features that we expect from 
it. The Lam- Tung relation, which is an exact property of the full leading-order contributions to the structure 
functions, independent of the coordinate frame, is not affected by resummation, even though the individual 
terms entering in it are affected. Furthermore, in the CS frame, the ratios A — {Wt — Wl)/{Wt + Wl) and 
V — 2Waa/(Wt + Wl) will not receive large corrections from resummation at small Qt- In particular, when 
restricting to the annihilation process, resummation has no effect on the cos 20 asymmetry. This observation 
is especially relevant for the Drell-Yan process in irp or pp scattering. The ratio fi = Wa/{Wt + Wl), on the 
other hand, will be subject to considerable resummation effects, due to the lack of the leading-logarithmic terms 
in its numerator. In the GJ frame all ratios A, /i, v will be affected by NLL resummation effects, which are not 
necessarily small. Similar conclusions were drawn for the ratios Ai. 

We finally emphasize again that we hope that our study will provide motivation for a development of the full 
NLL resummation for the structure functions Wl^a,aa- Given the experimental information available on the 
Drell-Yan process, on SIDIS, and on e+e~ annihilation, this would also have clear phenomenological relevance. 
Of further interest would be studies of the angular distributions integrated over Qt \^- Here, "threshold- type" 
logarithms may emerge (44> A5\ , whose structure and resummation in azimuthal distributions have so far not been 
investigated. 
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